Pseudopotential method for higher partial wave scattering 
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We present a zero-range pseudopotential applicable for all partial wave interactions between neu- 
tral atoms. For p- and d-waves we derive effective pseudopotentials, which are useful for problems 
involving anisotropic external potentials. Finally, we consider two nontrivial applications of the p- 
wave pseudopotential: we solve analytically the problem of two interacting spin-polarized fermions 
confined in a harmonic trap, and analyze the scattering of p-wave interacting particles in a quasi- 
two-dimensional system. 
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PACS numbers: 34.50.-s, 32.80.Pj 



Modeling of two-body interactions is the basic step in 
development of theories of many-body systems. In the 
ultracold regime, atomic collisions are dominated by s- 
wave scattering, and interactions can be accurately mod- 
eled by the Fermi- Huang pseudopotential 1, 2]. The 
situation changes, however, in the presence of scatter- 
ing resonances, which can strongly enhance the contri- 
bution from higher partial waves. Such possibility has 
been demonstrated in recent experiments by employing 
Feshbach resonances to tune the interactions of identical 
fermions in p-wave 0, 0]. In this context, the develop- 
ment of pseudopotentials for higher partial wave scatter- 
ing is of crucial importance for the theoretical description 
of ultracold gases with I ^ interactions. 

There are several approaches in the literature to de- 
rive pseudopotential valid for all partial waves. The first 
derivation comes from Huang and Yang 0. Their pseu- 
dopotential, however, is incorrect with respect to I > 
waves, as recently shown inj5|]. Several alternatives have 
been proposed |j 0, 0,0, EJ, which have specific limi- 
tations. For instance, [5j entails no regularization and is 
only applicable to mean- field theories; |6| requires knowl- 
edge of the wave function in the inner region of the shell 
potential, and taking the limit of shell radius going to 
zero in the final step of the calculation. 

In this Letter we address the problem of interactions 
in all partial waves. We correct the original derivation 
of Huang and Yang and obtain a comparatively simple 
pseudopotential. Next, we derive explicit pseudopoten- 
tial forms for p- and d- wave interactions, that are are 
very convenient in calculations involving anisotropic ex- 
ternal potentials. We illustrate this by solving analyti- 
cally the problem of two identical fermions confined in 
an anisotropic harmonic trap. We finally turn to interac- 
tions of atoms in low dimensional systems. We apply our 
p-wave pseudopotential to analyze the p-wave scattering 
in quasi-two-dimensional (Q2D) system, and show the 
occurrence of confinement-induced resonances analogous 
to s-wave scattering 

nana, 

and p-wave scattering in 
Q1D [12|. Our analysis is of direct interest for studies of 
controlled interactions between tightly confined fermionic 
atoms, relevant for applications to quantum information 
processing. 



First, we derive the pseudopotential for interactions in 
all partial waves. We start from the Schrodinger equation 
for the relative motion 



^(A + fc 2 )*(r) = V(r)tf(r), 



(1) 



where k 2 — 2[iE/h 2 , and /i denotes the reduced mass. 
We assume that the potential V(r) is central and has a 
finite range. Outside the range of the potential, the wave 
function exhibits the following asymptotic behavior 



(2) 
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where Yi m (0, (f>) are spherical harmonics, the radial wave 
functions Ri = Ci m [ji(kr) — tanSini(kr)] are linear com- 
binations of spherical Bessel and Neumann functions 
ji(r) and n/(r), C/ m are coefficients that depend on the 
boundary condition for r — > oo, and the phase shifts 
Si are determined by the potential V(r). Our goal is 
to replace V(r) by contact potential V ps (r), which acts 
only at r = and gives the same asymptotic func- 
tion v I / a(r) as the real potential V(r). Following Huang 
and Yang 0, we determine the pseudopotential from 
V ps (r)* a (r) = h 2 /(2^)(A + fc 2 )* Q (r). Since j,(kr) is 
regular at r = 0, the only contribution to the pseudopo- 
tential comes from the behavior of ni(kr) for small r: 
ni{x) (21- Hence 



V ps (r)y a (r) =— J2 (2^-l)!!C im tan ( 5 i y ir 

" 1=0 rn=-l 
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To calculate the action of the radial part of the Laplace 
operator on the function \/r l+1 at r = 0, one has to re- 
sort to the theory of generalized functions. By applying 
the Hadamard finite part regularization of the singular 
function l/r fc [3 El, one can prove the following iden- 
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where S^ n '(r) denotes the n-th derivative of the delta 
function. In comparison, Huang and Yang calculate (0 
by mapping on S(r)/r l , which leads them to the in- 
correct result. In the final step of the derivation we 
express the coefficients C; TO in terms of the regular- 
ization operator d 2l+1 r t+1 : C lm = (21 + l)!!/[fc'(2Z + 
iy.][df l+1 T l+1 Ri(r)] r=0 . In this way we obtain the fol- 
lowing form of the pseudopotential: 



Vps (rMr) =J- £ [(2/ + 1)!!]2 W ' 6(l){r) 



2[i 
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dr ,2l+i r 
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dQ,'Pi(n-n')^(r') 



r'=0 

(5) 



where the angular integral over O' acts as a projection 
operator on a state with a given quantum number I. 
Here, Pi(x) is the Legendre polynomial, n = r/r, and 
n' = r'/r'. We note that in the calculation of the matrix 
elements of V ps (r) the differentiation of the delta func- 
tion is equivalent to differentiation of the function that 
acts on the l.h.s. of the pseudopotential, with a proper 
change of sign. Moreover, for functions behaving like 
r for small r, one can substitute S^(r) = (-l)H\5(r)/r l 
and S(r) — 4irr 2 S(r) which shows that the s-wave compo- 
nent of 10 is equivalent to Fermi-Huang pseudopotential: 
— 2ttH 2 tan 8q /(fik) 5 (r), whereas components with / > 
differ from the pseudopotential of Huang and Yang by 
a prefactor. 

Now we derive an alternative form of the pseudopo- 
tential, with projection operators expressed in terms of 
the differential operators. Such representation is partic- 
ularly useful for problems involving anisotropic external 
potentials, since in this case the wave function, contain- 
ing several components of the angular momentum, has 
to be projected on a given partial wave interaction. Let 
us first focus on p-wave interactions. Expressing radial 
derivatives in terms of the gradient operator d r — n- V, 
where V = (d x , d y ,d z ), and performing the angular inte- 
gration in the projection operator, we obtain the follow- 
ing form of the pseudopotential for p-wave interactions 



V P (r) 



nh 2 a p (k) 3 
A* 



-a 3 

V 6(r) V r^i 



(6) 



where a p is the p-wave scattering length: a p (k) 3 — 

— tanSi(k)/k 3 , and the symbol V (V) denotes the gra- 
dient operator that acts to the left (right) of the pseu- 
dopotential. To derive 10, in the operator acting to the 
r.h.s. of the pseudopotential, we introduced an addi- 
tional multiplication by r: d 3 r 2 — > d r rd 3 r 2 = n-'Vrd^r 2 , 



which is equivalent for terms of the order of r, giving 
the contribution to the matrix elements. In this way we 
preserve the form of the regularization operator, which is 
crucial for an exact treatment of the interacting atoms. 
We note that the pseudopotential for p-wave interactions 
in the form containing a scalar product of gradients has 
been derived for the first time by Omont Q, however, 
in [3 it contained an incorrect prefactor, and did not 
include the regularization operator. We stress that the 
presence of the full regularization operator d 3 r 2 on the 
r.h.s. of the pseudopotential is important, since for 
anisotropic external potentials the exact wave functions 
obtained in the pseudopotential method contains terms 
behaving like Xi/r for small r. Such terms, for instance, 
are not removed by the operator V<9 2 r 2 

A similar procedure can be repeated for d-wave inter- 
actions. In this case we obtain the following pseudopo- 
tential 



V d (r) 



3d D v 

i,j,k,l 



u(dl x )5(v)(d 2 Xk 



r r , 

dr 5 ' 



(7) 



where D ijk i = S ik 5ji - f\S tj S k h 9d = Trh 2 a d (k) 5 /(8/z) 
and a,d(k) is the d-wave scattering length: ad(k) 5 = 
— t&nS2(k)/k 5 . As can be easily verified, the pseudopo- 
tential J7J) does not give any contribution for functions 
exhibiting s- and p-wave symmetries, which results from 
the implicit projection on d-wave contained in (0. 

Two interacting spin-polarized fermions in a harmonic 
trap. — We start from the integral form of the Schrodinger 
equation for the relative motion of atoms 



*(r)= d 3 r' G(r,r')V p (r')V(r'), 



(8) 



where G(r, r') = (r\(E — H) 1 |r') is the single-particle 
Green function, and H is the Hamiltonian including the 
external potential. For a harmonic potential G(r,r') can 
be represented by the following integral 



G(r,r') = - 



(VxVyVz) 1/2 

(2nf/ 2 d 3 hw J 

x k x 



n 



cxp \r/ k - 



dte Et/(hu) 



■h( x l+ x 'k) cosh ( t Vk)'' 



d' 2 sinhfir/fc ) 



sinh^fct) 1 / 2 



(9) 



which is convergent for energies smaller than the energy 



of zero-point oscillations E = h(u x 



-lu z )/2. FotE> 



Eq the Green function can be determined by analytic 
continuation of 0. Here, d = \Jh/(puj), r/k — ^ 
denotes some reference frequency, and the product in (0 
runs over k = x,y, z. 

Inserting the pseudopotential into (0, we 
obtain a set of linear equations J2i^ki(E)ci — 
Cfe, where the coefficients c t = \d Xl rdlr 2 ^>(r)\ r 
and the matrix elements are given by A k i (E) = 
ir(a p /d) 3 {d Xk rd 3 r 2 d x >G(r,r')] r=r/= o- The energy levels 
are determined from the secular equation of the matrix 



3 



A(E) — I. Below we discuss the case of axially symmetric 
traps: lu x = u> y = u>± = r/u> z . For pancake-shape traps 
with the anisotropy coefficient r\ = 1/n where n is some 
integer, the equations determining the energy spectrum 
have the relatively simple form 
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n—1 y( fc+1/2 



n p/ fc+l/2 
fc=0 1 ^ n 



?71 = 0, 



(10) 



,2 



fc+z+1 



m = ±l, (fl) 



where 5 = E/(Hui z ), d z — y/H/(fiu> z ) and to is the projec- 
tion of the angular momentum on the z-axis. As it can be 
easily verified, for the case of spherically symmetric trap 
= 1), Eqs. (|Tn|) - ((TT|l coincide with the result of Refs. 
Q . We note that in the presence of an anisotropic the 
degeneracy between states with m — and m = ±1 is 
lifted. We have derived similar closed analytic formulas 
for cigar-shape traps with 77 = n however, we do 
not present them in the present letter because of their 
rather complicated form. In general, when 77 7^ 1/n and 
?/ 7^ 71, the energy levels can be determined from an im- 
plicit equation involving the integral representation @ 

m 

Figs. Ola) and ^b) show the energy levels in har- 
monic traps with 77 = 1/10 and 77 = 10, for different 
values of the quantum number to. We observe that for 
small and positive values of the scattering volume, the 
system does not have eigenstates with energies smaller 
than the lowest eigenenergy for noninteracting atoms, 
This particular behavior results from the properties of 
the p-wave pseudopotential. In this regime of energies, 
however, it is necessary to include in the calculation the 
energy-dependence of the pseudopotential. To this end 
it is sufficient to insert in the pseudopotential @) the en- 
ergy dependence of the phase-shift Si(k), and calculate 
the energy spectrum in a self-consistent way. This ex- 
tends the applicability of the pseudopotential method to 
the large scattering lengths and tight confining potentials 
and allows to describe the entire molecular spectrum of 
the realistic potential 0, 0] . 

To illustrate this procedure we calculate the energy 
spectrum for a square-well model interaction, assum- 
ing for simplicity spherically symmetric harmonic trap. 
Fig. ^c) presents the energy levels of p-wave interact- 
ing particles as a function of the square-well depth U 
calculated for the square-well radius Rq = 0.05d. The 
depth U is varied close to the resonance scattering for 
p-wave, related to the appearance of a bound state 
at U = — 1974fiix>. The exact energy levels are com- 
pared with predictions of the pseudopotential method 
given by Eq. (|1(J|> . and with the energy spectrum cal- 
culated in a self-consistent way from Eq. i|10[l with a 3 
replaced by a p (k) 3 = — t&nSi(k)/k 3 . For comparison, 
Fig. ^d) shows the energy levels of s-wave interacting 
atoms, given in the pseudo pot ential approximation by 
d/a s = 2T(§ - f )/r(i - |) 18], where £ = E/{hu) and 
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FIG. 1: (a)-(b): Energy spectrum of the two fermions inter- 
acting in p-wave as a function of the scattering volume a^/d 3 , 
with d — -\/h/fiuj z . The two particles are confined in an axi- 
ally symmetric harmonic trap with rj = uj±/uj z =0.1 (a), and 
77 = 10 (b). (c)-(d) Energy spectrum of two atoms interacting 
via a square-well model potential as a function of the well 
depth U. The atoms are confined in an isotropic harmonic 
potential and interact in p-wave (c) and s-wave (d). 



a s is the s-wave scattering length. 

One observes that self-consistent calculation with 
energy-dependent pseudopotential (EDP) provides very 
accurate results for the energy spectrum. On the other 
hand, the ordinary pseudopotential method fails for large 
scattering volumes, and for energies where the bound 
state of the square well potential appears. Moreover, its 
range of applicability decreases for higher energy levels. 
For comparison, the ordinary s-wave pseudopotential is 
incorrect only with respect to deep bound-states. 

This behavior can be explained by analyzing the effec- 
tive range expansion: fc 3 cot <5i(fc) = — 1/a 3 — k 2 / (2R*) + 
0(k A ), where R* can be interpreted as effective range for 
p-wave (for square- well R* — i?o/3). The second term in 
the expansion can be neglected when fe|a p | <C (kR*) 1 ^ 3 , 
which combined with the condition for the applicabil- 
ity of the pseudopotential: kR* <g; 1, gives fc|a p | <C 1. 
Therefore, the regime where one can apply the ordinary, 
energy-independent pseudopotential for p-waves is quite 
narrow. This can be compared to the s-wave, where the 
analogous condition takes form k\a\ <C (kR) -1 , where R 
is the effective range for s-wave scattering and (fci?) -1 is 
large. 

Scattering of spin-polarized fermions in Q2D 
systems. — The scattering solution can be found 
from the Lippmann-Schwinger equation: 

*(r) = * (r) + Jd 3 r' G + (t,t')V p (v')*(v'), (12) 

where 5'o( r ) represent the wave function of the incoming 
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FIG. 2: Two-dimensional differential cross section in forward 
direction |/(0)| 2 in the scattering of two fermions in a quasi- 
2D system as a function of p-wave scattering volume a^, for 
different relative energies E of the particles. The inset shows 
the value of the critical volume V c , at which the system ex- 
hibits the confinement-induced resonance. 



particle and G+(r,r') = (r|(£ , -H + ie)" 1 |r / ). The Green 
function G+(r, r') describing the propagation of outgoing 
waves, can be determined from Eq.@ by taking the limit 
luj_ — > and performing the analytic continuation for 
energies E > E . 

In 2D the scattered wave in the asymptotic regime 
(p — > oo) is described by /(</>)e lfep / y/p with the scattering 
amplitude /(</>) = {2mk)' 1 / 2 E~=_^ e^e* 2 *™ - 1), 
where 5 m are scattering phase-shifts [19j. In the regime 
of energies hw z /2 < E < 3hui z /2 the motion in the z 
direction is frozen and sufficiently far from the scattering 
center it is described by the ground-state wave function 
of the harmonic oscillator: tpo(z) — e~ z /( 2< **) /V 1 / 4 . By 
solving Eq. $Tty with the pseudopotential ©, we obtain 
the scattering solution, which for p-wave interactions con- 
tains only m — ±1 scattering waves 



p— »oo 



■00 (2) ^ e 



3 ikp 



4cos< 



1 + i cot <5i y/2mkp 
The scattering phase shift Si is given by 



(13) 



cot Si = — 



(14) 



where W(x) is a function taking values of the order of 
one [52, and £ = E/ {hw z ). When cot Si = 0, the system 
exhibits confinement-induced resonance, which occurs at 
the scattering volume V c (£) = y/^d 3 z /W{£/2 - 1/4). At 
low energies {£ — ► i) function W(a;) has a well defined 
limit: W(0) w 0.325, and V c (0) « 5.4^. This qual- 
itatively differs from s-wave interactions, where cot <Jq 
exhibits logarithmic behavior for small k 01 an( i the 
resonance condition for k — > depends on the relative 
kinetic energy of the scattered particles. Fig. [21 shows 
the dependence of the differential cross section \f((j>)\ 2 at 
4> = on the scattering volume a 3 for different values 
of the energy. We observe that at low energies of the 
scattered particles the curve is strongly peaked around 
V c . For higher energies the resonance is broader and fi- 
nally for £ sa 0.525 disappears. In this regime of energies 
the resonance is present for positive values of the scatter- 
ing volume, which can be observed in the inset of Fig. [21 
presenting V c {£). 

In summary, we presented the zero-range pseudopo- 
tential applicable for all partial wave interactions. For p- 
and d-waves we derived an alternative representation of 
the pseudopotential, in which the projection on spheri- 
cal harmonics is replaced by an appropriate differential 
operator. The p-wave pseudopotential has been applied 
to calculate analytically the spectrum of two interacting 
fermions in a harmonic trap, and to study the scattering 
of identical fermions in a quasi-two-dimensional system. 

After completing this paper we learned of recent work 
of Derevianko |20j in which the pseudopotential equiva- 
lent to JSJ is derived. 
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